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O^l . We investigate the distance-redshift relation in a realistic inhomogeneous universe where 

the mass distribution is described by the mass function of Sheth and Tormen. It is found 
that the derived distance deviates systematically from the standard distance up to 10% 
depending on the choice of the lowest halo mass in which baryonic matter condensed to form 
luminous object such as galaxies. Remarkably the derived distance is well approximated by 

■ the Dyer-Roeder distance if we choose the dumpiness parameter a calculated by our model. 
We also discuss the effect of inhomogeneities in the determination of dark energy parameter 
in the supernovae observation, and find that this effect must be taken into account for the 

■ future high redshift supernovae observation. 



§1. Introduction 



It is needless to say that the distance-redshift relation has a fundamental im- 
portance in the observational cosmology. The distance formula is usually derived by 
assuming that the universe is homogeneous and isotropic. This assumption has been 
supported by the observations of large scale galaxy survey and cosmic microwave 
background radiation over the scale of about lOOMpc. The derived formula is called 
the standard distance or the filled beam distance. However the light rays feel not 
the averaged geometry but local inhomogeneous geometry, and thus it is reasonable 
, to expect that the local inhomogeneous matter distribution does have some effect 

on the distance-redshift relation. There have been various attempts to derive the 
distance-redshift relation taking the effect of dumpiness of matter distribution into 
account. 1 ^ 6 ) Most well known is the so-called Dyer-Roeder (DR) distance which 
is derived under the assumption that the light rays propagate the vacuum region 
along their entire paths. This is a reasonable assumption in a sense since we do not 
observe the distant light sources inside galaxy this side of the source. In fact there 
are some studies supporting the validity of DR distance in a simplified situation 
where dark matter distribution is not taken into account. However it is now widely 
accepted that dark matter is an important constitution of the universe and plays an 
essential role in the structure formation. Thus the correct distance-redshift relation 
must take the distribution of dark matter into account. The purpose of this paper 
is to construct such a distance-redshift relation. 

In this paper we describe the dark matter distribution in terms of the mass 
function proposed by Sheth and Tormen. Moreover we assume that the light rays 
from the distant light sources cannot propagate through the mass concentration 
above a minimum mass. This is supposed to mean that the the baryonic matter 
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collapses to form luminous objects such as galaxies in these mass concentrations 
and prevent the distant light. Thus the light rays propagate statistically underdense 
regions and it is expected that this will cause systematic deviation from the standard 
distance-redshift relation. 

The structure of this paper is as follows. We first describe our model of inho- 
mogeneities and explain our assumption in the light propagation. In Section 3 we 
then apply the above model to the distance-redshift relation given by Futamase and 
Sasaki which applies realistic inhomogeneous universe. We will find that the relation 
systematically predicts larger distance than the standard distance depending on the 
lowest mass introduced in Section 2. We will find in Section 4 that the so-called 
Dyer-Roeder distance gives a very good approximation to the derived distance if we 
use the dumpiness parameter a calculated by our model of inhomogeneities. We 
apply our result to the determination of the parameter of equation of state in the 
dark energy in Section 5. Finally Section 6 will be devoted to some discussion. 

For fiducial cosmology, we adopt in this paper a totally flat universe model with 
h = 0.7, ^ de = 0.72, f2 m = 1 - Q de = 0.28, Q h = 0.046, w = -l,w a = 0,cr 8 = 
0.82, n s = 0.96. 

§2. Model of mass distribution and the light propagation 

In a real universe, matter distribution is very clumpy and we model this clumpy 
universe by halo model. 7 ) The comoving number density of collapsed object (halo) 
which has mass m in the mean comoving mass density p m at redshift z, n(m,z), is 
described by 

m 2 n(m, z) dm ,, .dv , 

±-L-L— = V f(y)— 2-1 

p m m V 

The height of density peak, is, is defined as 

v = 2 C V \ 2-2 
a z (m, z) 

where S c (z) ~ 1. 686 {2 m 0055 (z) is the overdensity of a spherical region that collapses 

Sit 

(z) is the ratio of matter density to critical density at z. a(m, z) is 
the rms value of the density contrast when smoothed with a tophat filter of scale 
R = (3m/47rp m ) 1 / 3 , extrapolated using linear theory to z : 

a\m,z) = Jf k3P 2f ,Z) \W(kR)\ 2 (2-3) 

where P lm (fc, z) is the linear matter power spectrum at z and W{x) = (3/x 3 )(sina; — 
x cos x) is tophat filter in fourier space. To predict the halo number density, we use 
the Sheth and Tormen's mass function 9 ) 

■*M = A ( 1+ ;»)U) V (24) 
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with v' = 0.75is,p = 0.3 and normalization constant A is chosen to satisfy the 
constraint f f(u)du = 1. When we calculate the linear matter power spectrum, we 
use the BBKS transfer function, 10 ) the shape parameter from Sugiyama 11 ) and the 
growth suppression rate, g(z) = Dy m (z) / a(z) , which obey the following differential 
equation 

+ I 5 " 3 ^^ dc ^dh7^ + ^ - ^)}f2 de (z)g(z) = (2-5) 

where a(z) = 1/(1 + z) is the scale factor, f2 de (z) is the ratio of dark energy density 
to critical density at z and w(z) is the equation of state of dark energy. Hereafter 
we parametrize w(z) = w o + w a (l — 1/(1 + z)). 

In the standard structure formation theory, halo which has sufficient mass can 
attract baryonic matter and forms luminous objects such as galaxies. So if we chose 
this mass scale appropriately, we can predict the number of halo which contains 
luminous objects and mass of these halos. If these halos exit between observer and 
sources, light rays from sources have the possibility to be prevented by these halos 
or extincted by dust associated with these halos. In this paper, we consider this 
possibility to its maximum. For very large halo, like a cluster of galaxies with very 
large M/L ~ 200, this assumption may not be appropriate, because of large fraction 
of the 'transparent' matter. But the effect of these large halo is very small in the 
distance as we will see later. 

Based on the above assumption, we can define the effective comoving matter den- 
sity (or smooth distributed matter density) which light right ray can passes through 

Pes{z)=p m - dmmn(m,z) (2-6) 

JM min 

M m i n is a minimum halo mass in which baryonic matter condensed to form 
luminous objects and is our parameter. Figure 1 shows the ratio of effective comoving 
mass density to total comoving mass density with various choice of the minimum 
mass. One can see that very large halos (M m i n = 1O 14 M ) have little effect on 
the effective density. This is because they have large masses but is very small. On 
the other hand, halos with mass range 10 n M Q < M < 1O 13 M contributes most 
effectively because m x n(m, z) is the largest around 1O 12 M at z = 0. As redshift 
gets larger, an inhomogeneity becomes small and p c s approaches to p m at z ~ 0(1). 
In these calculation, we chose highest masses infinity since the contribution by masses 
much larger than 1O 14 M may be totally neglected. 

§3. Distance-redshift relation 

In a homogeneous flat FRW universe, the angular diameter distance is defined 
as 12 ) 



(3-1) 
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M > 10" M 



0.01 0.1 1 10 

z 

Fig. 1. The ratio of effective comoving matter density to total comoving matter density as 
a function of redshift for various mass range: lowest mass is selected at 10 n M© (solid), 
1O 12 M (dashed), 10 13 M Q (dot-dashed) and 10 14 M Q (dotted). 



x( z ) is the comoving distance in a homogeneous universe, which is an integral over 
redshift of inverse Hubble parameter H(z) 

Using the optical scalar equations 13 -* and the perturbation theory, Futamase 
and Sasaki showed that in a realistic inhomogeneous universe, the angular diameter 
distance is modified as 14 )~ 16 ) 



x / M Sd F A RW (z,h) iA 1 



1 



(vlhi - v l hi) 



aH 

/ ax (1 + z)5 m (z ,n) (3-3) 

Jo Xs 



where h is the source direction, Xs is the source comoving distance, v s , v is the source 
and observer peculiar velocity respectively and 5 m is the matter overdensity. The first 
line is so called the doppler term because of local peculiar velocity, which changes the 
redshift of sources relative to observer and the solid angle of observer. The second 
line is so called the lensing term because of inhomogeneity of line-of-sight matter 
distribution (more detailed physical description is described in Ref. 17)). Then, we 
define 

df(z) = dl RW (z)(l + (5 dA (z,n)) b ) (3-4) 

(. . . ) means ensemble average and subscript '6' denotes that we can only average 
using observed sources. Of course, if we can average all perturbed density field, this 
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distance return to c^ RW . But under our assumption discussed in §2, average of 5 B 
doesn't vanish 



(K{z,n)) b = M _ i 

Pm 



(3-5) 



This means that light ray which can arrive us feels limited overdensity and tends to 
pass through underdense region, so doesn't equal to (i^ RW . From Fig.l, p e s/p m 
is always smaller than 1, so we expect that <i^ s is always larger than g^ rw . 

Figure 2 shows the fractional differences of and d^ RW as a function of redshift 
for various mass range. As expected above, d^ is always larger than <i^ RW . The 
lensing weight function, (xs ~ x)x/Xs, makes this effect continue to grow by higher 
redshift compared to (5 m ) b . Very large halos (M > 10 14 M©) change distance a little 
(< 0.1%) and we can ignore. But middle large halos change distance largely up 
to 10% and we can't ignore at high redshift apparently. At z ~ 1, this effect look 
mildly small, but in searching for dark energy property precisely this change may 
make some bias to dark energy parameter, we will discuss in §5. 




Fig. 2. The fractional differences on angular diameter distance between rf^ an< i o!a RW as a function 
of redshift, which is defined as (dI s /dI RW - 1). They are defined in text respectively. The mean 
of curves is the same as Fig. 1. 



§4. Dyer-Roeder distance 

There have been many discussions on a possible effect of local inhomogeneities 
in real universe on the distance-redshift relation starting from the discussion by 
Zel'dovich. 18 ) For example, Dyer and Roeder derived the distance formula on the 
assumption that a certain fraction, the constant dumpiness parameter a, of all 
matter is distributed uniformly whereas the rest is clumped into galaxies, and that 
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light rays travel well away from all clumps of matter, feeling only the effect of the 
fraction a of all matter. 19 ) In Dyer-Roeder formalism, angular diameter distance, 
d ] ^ R (z), obeys the following differential equation 



(1 + z fF^- + (1 + z)Q d ^- + HdT 

dz z dz 

which satisfies the boundary conditions: 

dg R (0) = 



dd» R 



dz 



2=0 



1 

#0 



where 
T=H 2 {z) 

G = 

_ 3a(z 



(l + z)H(z)^- + 2H 2 {z) 



3 | 3(w(z) + 1) 



• cxp 



! f if>'> +1 > 



(4-1) 

(4-2a) 
(4-2b) 

(4-3a) 
(4-3b) 

(4-3c) 



For our fiducial model, w(z) = — 1, the second term of TL vanishes. Of course, if 
a = 1, H becomes (1 + z)H(z)dH/dz and return to the FRW distance. In our 
inhomogeneous model, the dumpiness parameter a is equal to p e s/f>m as plotted in 
Fig 1. 

Figure 3 shows the fractional differences of (i^ R and c^ RW as a function of 
redshift for various mass range. For a < 1, light beam is defocused and c^ R is 
allways larger than d^ RW . The figure shows that it is remarkably close to with 
deviation < 0.3%. 



§5. Parameter biases for future SNIa survey 



In this section, we consider how this difference on distance affects the estimation 
of cosmological parameters based on the observational probe, especially Type la 
supernova (SNIa) survey. For SNIa survey, observable quantity is the apparent 
magnitude m(z) which relates the luminosity distance 

m(z) = 5\og 10 d h {z) + M = ^-^\nd L (z) + M (5-1) 

where M is the magnitude zero-point. Rescaling c?l, we can subtract Hubble con- 
stant Hq uncertainty by the redefinition of M. So our interest parameter is 0{ = 
{f?de, wq, w a , M}. Recalling the relation between the luminosity distance and the 
angular diameter distance, di, = (1 + z) 2 dA, we can write difference on apparent 
magnitude in our model 



^{z) = -^-\n{l + {5 dA (z,n)) b ) (5-2) 
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Fig. 3. The fractional differences on angular diameter distance between d A and d A as a function 
of redshift, which is defined as (rfA R /^A RW — !)• They are defined in text respectively. The mean 
of curves is the same as Fig. 1. 



SNIa survey has some uncertainty which comes from not only astrophysical, but 
also cosmological ones. 17 )' 20 ) In this paper, we only consider the intrinsic dispersion 
of magnitude a m . If the noise covariance consists of only <r m , bias to parameter 
vector 6i is described as 21 ) 

80i = F-\ 3 \ f dzN(z)Sm(z)^ (5-3) 
where Fy is the Fisher matrix of fiducial survey : 

and N(z) is the number of SNIa at z. We consider SNAP like survey and use survey 
parameters written in Kim et al. 21 ) Substituting (5.4) into (5.3), <r m is cancelled 
out. So we don't have to chose a m in bias calculation and our result is not affected 
by cr m choosing (If we consider other uncertainty, this does not hold anymore). 

Firstly, we consider bias to parameter vector 6i = {^de^Oi^} (constant dark 
energy model) to understand physical meaning. Table 1 shows biases to parame- 
ter estimation for various mass ranges with and without prior a(l7dc) = 0.03 which 
is anticipated from other cosmological probes. Intuitively, to make distance larger, 
i?de becomes larger or w$ becomes smaller. Table 1 makes this expectation apparent. 
Without prior, i?dc becomes larger, but wq becomes larger because of degeneration. 
With prior, fi& c becomes larger slightly because of prior limitation and wq becomes 
smaller for compensate larger distance. It is found that the effect of local inhomo- 
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geneities affects the parameter wo of the order of 0.1, thus this effect must be taken 
into account in the future supernovae survey. 

Table I. The biases to parameter vector 0i — {i?dc, wo, M} for various mass ranges with and without 
prior a(Q dc ) = 0.03. 



mass range Sf2 c \ c (w.o/w prior) Swq(w.o/w prior) 



M > 10 n M o 


7.4 x 10" 


78 


8 x 10" 


-3 


1.6 x 10" 1 / - 


1.8 x 10 _i 


M > 1O 12 M 


4.4 x 10" 


2 /5 


2 x 10" 


-3 


8.2 x 10" 2 / - 


1.2 x 10" 1 


M > 1O 13 M 


1.5 x 10" 


Vi 


8 x 10" 


-3 


2.1 x 10" 2 / - 


4.9 x 10" 2 


M > 1O 14 M 


1.3 x 10" 


Vi 


5 x 10" 


-4 


-4.3 x 10~ 4 / - 


- 6.4 x 10" 3 


lO n M < M < 10 13 M Q 


5.9 x 10" 


2 /7 


1 x 10" 


-3 


1.4 x 10" 1 / - 


1.4 x 10" 1 



Next, we consider bias to parameter vector = {Q^ e ,wo,w a , M} (time variable 
dark energy model). Table 2 shows biases to parameter estimation for various mass 
ranges with and without prior cr(,!?de) = 0.03. Similar to above, to make distance 
larger, l?de becomes larger or wq becomes smaller or w a becomes smaller. Without 
prior, <5J?de and 5wo is always positive but 8w a is positive at low mass model and 
negative at high mass model. With prior, 5f2^e an d Swq are always positive and 5w a 
is always negative. We can understand behavior of <5i?de an d ^ w a similarly to above. 
But in spite of mass range and prior, 5wo is always positive. 

Table II. The biases to parameter vector 0i = {fide, wo, w a , M} for various mass ranges with and 
without prior a(Odc) = 0.03. 



mass range 5J?dc(w.o/w prior) 5wo(w.o/w prior) <5w a (w.o/w prior) 

M > 10 n M Q 1.0 x 10-73.4 x 10 -4 1.1 x 10 _1 /2.6 x 10" 1 5.5 x 10"7 - 1.4 

M > 1O 12 M 5.3 x 10" 2 /1.8 x 10" 4 6.7 x 10" 2 /1.4 x 10 _1 1.8 x 10"7 - 8.3 x 10" 1 

M > 1O 13 M 1.4 x 10" 2 /4.5 x 10 -5 2.3 x 10~ 2 /4.3 x 10" 2 -3.3 x 10~ 2 / - 3-0 x 10" 1 

M > 1O 14 M 3.9 x 10 _4 /1.3 x 10 -6 1.0 x 10- 3 /1.6 x 10~ 3 -1.8 x 10" 2 / - 2.6 x 10 -2 

10 11 Mq < M < 1O 13 M 6.9 x 10 _2 /2.9 x 10 -4 8.7 x 10 _2 /2.2 x 10" 1 3.1 x lO -1 / - 1.1 



We also consider influence on the current SNIa data, especially 'UNION' data. 22 ) 
Doing similar calculation above (we use each a m available from above URL), we 
find that bias to parameter vector 6i = j,!?^, i«cb -^} with prior o"(i?de) = 0-03 in 
mass model 10 n M o < M < 1O 13 M is 5Q Ae = 1.1 x 10~ 5 and 5w = -4.1 x 10~ 2 . 
Thus this effect is negligible compared to statistical error in the current SNIa data 
because of low redshift SNIa. 



§6. Conclusion and discussion 



In this paper we investigate the distance-redshift relation in a realistic inhomo- 
geneous universe where the mass distribution is described by the mass function of 
Sheth and Tormen. It is found that the derived distance deviates systematically 
from the standard distance up to 10% depending on the choice of the lowest mass. 
Remarkably the derived distance is well approximated by the Dyer-Roeder distance 



SNIa data is available from 



http://supernova.lbl.gov/Union/ 
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if we choose the dumpiness parameter a calculated by our model. We also discuss 
the effect of inhomogeneities in the determination of dark energy parameter in the 
super novae observation. The effect behaves like dark energy and thus is relevant in 
dark energy survey. Although it does not change the result of current SNIa survey, 
it should be seriously taken into account future SNIa survey since it aims at deciding 
the dark energy parameter in a few % statistical error. 

In the present formalism we have not considered the dispersion of the distance- 
redshift relation. It will be very interesting to consider the dispersion for more 
realistic interpretation of the observed data. The dispersion itself has important 
meaning in cosmology because it carries the information of large scale structure as 
well as dark energy property. The dispersion depends on the spatial and temporal 
distribution of halos and we would like to consider this problem in future work. 

We have applied our distance model to SNIa survey in this paper. It is also 
interesting to find other applications such as the evolution of the luminosity function 
of quasar based on our distance model. 
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